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(1) Determine a solução do seguinte problema de valor inicial: 

• tet + t ía\ n 

y = 7T-- — ' y ° = õ 

2 + smy 2 

Resolução: Uma vez que 2 + siny nunca se anula a equação é equivalente à equação 
separável 

(2 + smy)y = te t + t 
Integrando de a t obtém-se 



ou seja 



2y - cos y - 2y(0) + cos y(0) = te 1 - é + -if - (0 - 1 + 0) 



2y - cos y - tê + e t - -t 2 = n + 1 



Uma vez que 2 + siny(O) = 2 + sin \ = 3 ^ 0, o teorema da função implícita garante que 
esta equação define implicitamente y em função de t numa vizinhança de to = 0. D 



(2) Determine a solução do seguinte problema de valor inicial: 
(tety - 2y)y = -ye*» - 1 , y(0) = 1 

Resolução: A equação pode escrever-se na forma 

M{t,y) + N{t,y)y = 
onde M(t, y) = ye ty + 1 e N(t, y) = te ty - 2y. Tem-se 

— = e ty + tye ty = — 
oy ot 

portanto a equação diferencial é exacta. Um potencial 4>(t, y) para (M, N) é uma solução 
do sistema 

í W = yety + 1 ^ ! 4>(t,y) = e ty +t + A(y) 
\% = te ty - 2y ^ \ 4>{t,y) = e ty - y 2 + B(t) 

Uma solução é <p(t, y) = e ty + t — y 2 . Conclui-se que a solução do problema de valor inicial 
verifica a equação 

e ty + t - y 2 = 0(0, 1) = 
Uma vez que 
d 

(*,»)=(0,1) 



dy (e- + t-y 2 ) 



= ( teíy -^)la.„wo.i)=- 2 ^ 



(*,»)=(o,i) 

o teorema da função implícita garante que esta equação define implicitamente y como uma 
função de t numa vizinhança de to = 0. D 
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(3) Determine a solução do seguinte problema de valor inicial: 

.1 ty 

arctan y + e -\ õ = , vil) = 

y 2 + 2y 2 yv ; 

Resolução: Sejam M(t,y) = arctan y + e* e N(t,y) = 2 , 1 3 . Tem-se 



2+2y 2 



dM 1,1 diV 



^ 



<9y 1 + y 2 ^ 2 + 2í/ 2 dí 

/ogo a equação não é exacta. 

Para que exista um factor de integração p, = p,(t) é necessário que a seguinte equação 
tenha solução: 

d ( t 2 \ d ( tp,{t) 



-(MíXarctany + e* 2 )) = ^ (j^ 



M(*)^— — õ = M*) , o 9 + /*'(*); 



1 + y 2 ^ w 2 + 2í/ 2 ^ w 2 + 2í/ 2 



2 + 2y 2 ^ w 2 + 2í/ 2 

fj/(t)t = fi{t) 

Assim pode-se tomar p(t) = t para factor de integração quando í ^ 0. Multiplicando a 
equação por t obtém-se a seguinte equação exacta, que é equivalente à inicial para t ^ 0: 

t 2 t ■ r. 

t arctan y + te -\ ^ -y = 

Um potencial para esta equação é uma solução do sistema 

í §£ = t arctan y + te f2 \ 0(í, y) = y arctan y + ^e* 2 + A(y) 

1 I = 2W 1 0(*,2/) = f arctan y + 5 (í) 

.Ass/m uma solução é dada por 

,1 s f2 l t 2 

9{t,V) = — arctan y+ -e 



Como 0(1, 0) = ^e, a solução do problema de valor inicial verifica 

t 2 1 ,2 1 

yarctany + -e = -e 

t 2 
e — e 



arctan y — 

í 2 



y(í) = tan 



e — e* 



í 2 y 

O intervalo de definição é o maior subintervalo do conjunto 

D = {tem \ { }: -I < e —^ < 1 } 

que contém í = 1. fsíe intervalo não se consegue achar explicitamente uma vez que não 
se consegue resolver analiticamente as equações 

t 2 
e — e 7T 

~ í 2- 2 

D 
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(4) Determine a solução geral da seguinte equação diferencial: 

y - 2e~ y = 2te~ y 

Resolução: Esta equação é separável: 

y-2e~ y = 2te~ y 

e y y = 2í + 2 

e y = t 1 + 2í + C com C el 

y(í) = ln(í 2 + 2t + C) comC el 

O intervalo de definição de uma solução é um intervalo máxima! contido em 

{teR:t 2 + 2t + C>0} 
Uma vez que 



t 2 + 2t + C = ^> t= ^ ^> í = -1± y/l-C 

conclui-se que para que uma solução tenha domínio não vazio, tem que ser C < 1 e nesse 
caso, os intervalos máximos de definição das soluções são 



] -00,-1- Vi- C[ ou } -1 + Vl-C,+oo[ 

D 

(5) Mostre que o problema de valor inicial 

y = y* , 3/(0) = o 

tem infinitas soluções. Porque é que isto não contradiz o teorema de Picard? 

Resolução: Começamos por observar que o problema de valor inicial tem a solução 
constante y(t) = 0. Para j/^Oa equação é equivalente a 

y~*y = i 



2 U s ~Vo ) =t-t 
2 2 2 

2/ f = g(*-*0) + tfÔ 

2/(í) 2 = Q(*-*o)+v! 



//(O = ± \/ ( ^(* - *o) + í/l) P ara ^(* " *o) + vi > 



L/ma vez que 



u(t) = ±\/( g(*-*o) + l/l) parat>t --yl 



d ( (2 2\ã\ /2 i" 

5 ± U<*-«.) + ri) =±V 5 « -«.) + »! 
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tem-se 

lim 2 ^(t) = 
í->t -fs/,f 

2 

e portanto estas soluções podem prolongar-se ao instante t = to — |y 3 como funções de 
classe C 1 . Em particular, para yç, = tem-se "soluções" (as aspas referem-se ao facto de 
não estarem definidas num intervalo aberto) 



2 



v(t) = ±dl-(t-to)\ para t> t 

que se podem "colar" à solução constante igual a para obter soluções da equação definidas 
para íel: 

y (t) = \ ± VuK*-V)? Para t > t 
{ se t < t 

Se to > estas funções são soluções do problema de valor inicial do enunciado 1 . Conclui-se 
que o problema de valor inicial tem infinitas soluções. 

Isto não contradiz o teorema de Picard porque a função f(t, y) = ys não é localmente 
Lipschitziana em nenhuma vizinhança U de um ponto (ío,0). Se f fosse localmente 
Lipschitziana existiria C tal que 

\f(t,y)-f(t,0)\ ^ 

: < C para (t, y) e U 

\y - °l 

Mas em qualquer tal vizinhança 

U-0 



sup ■ 



|y-o| 



— - 

= sup \y s = +oo 



D 



(6) Esboce o campo de direcções e trace os respectivos tipos de solução para as seguintes 
equações diferenciais: 
(a) y = siny 

Sugestão: Na alínea b), comece por achar as soluções da equação 
da forma y(t) = et onde c é um número real. 

Resolução: 

(a) Para cada c 6 IR, o conjunto dos pontos (t,y) e IR 2 onde o gráfico da solução y(t) 

dy 
dt 



tem declive -# = c, é determinado pela equação 



sin y = c 



Não é difícil provar que estas são íocías as soluções. 
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Portanto os declives possíveis têm que estar no intervalo [—1,1] eo campo de direcções 
é invariante mediante translações de 2tt na direcção do eixo dos yy. 



Casos especiais: 



c= -1 

r- — _I 
2 



c = 



1 



2 

c = 



siny = — 1 



y= -| + 2fc7r 



smy = -2 <S= 
y = f + 2fc7r ou y = 
y = kn 
y=^ + 2kn 



y = -| + 2fc7r ou y = -^ + 2fc7r 



5tt 
6 



2fc7T 



C= 1 

Esboço do campo de direcções 



*y 



Traçado dos tipos de solução: 





(b) Para cada c e M, o conjunto dos pontos (t, y) e 
tem declive -^ = c, é determinado pela equação 



onde o gráfico da solução y{t) 



y + zt 

y-3t 



= c . 



e a equação não está definida para y = 3í. Note-se no entanto que quando uma 
solução da equação tende para um ponto sobre esta recta o declive do seu gráfico 
tende para infinito. Podemos portanto pensar nestes pontos como aqueles em que a 
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recta tangente às soluções é vertical. Em geral, tem-se 

V + 2t = 
Í/-3Í 

<í=> y + 2t = cy — 3ct 

^> (l-c)y = -(3c+2)í 

Portanto os conjuntos de pontos onde o declive das soluções é constante igual a c 
são rectas que passam pela origem. 

É natural ver se algumas destas rectas são soluções (cf. sugestão). Para c ^ 3, y = et 
é uma solução da equação sse 

ct + 2t 



ct- 


-3í 






c + 2 = 


c 2 - 


-3c 


c 2 - 


4c 


-2 


= 



^^ c = 2±V6 
Casos especiais: 

C=l É=0 

c=-| y = 

c = y = -2t 

Esboço do campo de direcções: 

t y / y=3t 

/ / / 

V;// 



/ ' 

/ 
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Traçado dos tipos de solução: 




D 



(7) Para cada uma das seguintes matrizes A ache uma forma canónica de Jordan J e 
uma matriz de mudança de base S tal que A = SJS' 1 . 



3 
-2 2 



~ 3 
2 


í ~ 

2 


1 

2 

3 


5 
2 . 

1 " 


-2 


2 



(a) A = 
(b)A = 
(c)A = 
(d)A = 

(e)A = 



Resolução: 

(a) Os valores próprios de A são 



0" 




0-1 




1 




0-10" 


-1 -1 


1 






3- A 

-2 2- A 



= ^^ (3- A)(2- A) = 
portanto a matriz é diagonalizãvel e tem-se 

J = 
Os vectores próprios para A = 2 são os que verificam 

^^ a = 



2 
3 



A = 3 ou A = 2 



1 

-2 



a 
b 
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portanto uma 


base dos vectores 


próprios 


é constituída por 


v\ = 



1 




Os vectores próprios para A = 3 são os que verificam 




" 
-2 -1 




a 
b 


= 


" 



^^ b= -2a 



portanto uma base dos vectores próprios é constituída por, por exemplo, 

1 
V2= [-2 

A matriz de mudança de base é 

S 

(b) Os valores próprios de A são 

3 _ \ I 

= 

1 

- = 



~1 


1 




5 = 


1 -2 




) 

3 X 1 
2 " i , 2 

~2 2 ~ Á 


(|-A)(|-A) 


A 2 - 4A + 4 = 


(A-2) 2 = 


A = 


2 





Como não pode haver dois vectores próprios linearmente independentes (nesse caso 
o espaço próprio seria todo o M 2 e portanto A teria de ser igual a 21 o que não é o 
caso) conclui-se que 

' 2 1 



J = 



2 



Um vector próprio associado a X = 2 é uma solução da equação 



(A - 2I)vi 




= 






1 1 " 

1 1 

2 2 . 




a 
b 


= 


' " 



a = b 











portanto uma base dos vectores próprios é constituída por, por exemplo, 

1 



vi = 



1 



Um vector próprio generalizado associado ao vector próprio v\ é uma solução da 
equação 

{A - 2I)V 2 = V! 

_1 1 

I I 



2 2 

1 1 

-a + -6=1 

2 2 



Pode-se tomar, por exemplo, a = -1 e b = 1, isto é 

-1 



v 2 = 



1 
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Conclui-se que uma matriz de mudança de base S tal que A = SJS -1 é 

S = 



1 -1 

1 1 



(c) Os valores próprios de A são 



3- A 1 

-2 2- A 



= 



(3- A)(2- A) + 2 = 
A 2 - 5A + 8 = 



A = 



5± 



-7 



portanto a matriz é diagonalizável e tem-se 
J = 



5+iV7 n 

2 u 

5-iV7 
2 







Os vectores próprios para A = 5+ ^ são os que verificam 



l-iy/7 
2 



1 

íWZ 



a 
b 



b = 



i\/7 



portanto uma base dos vectores próprios para A = 5+l ^ 7 é constituída por 



v\ = 



1 

-lW7 



5-Vji 



Como A é uma matriz real, os vectores próprios para A = % s ^° os conjugados 

dos vectores próprios para A = 5+ / 7 % portanto uma base dos vectores próprios é 
constituída por 



V 2 = 


1 

9 






A matriz de mudança de base é 


S = 


1 1 

-l+iy/7 l+iV7 
L 2 2 




(d) Os valores próprios de A são as soluções de 




-A 






-A 1 


= 




1 -A 




^^ -A 3 -A = 


^^ A(A 2 + 1) = 


<í=> A = ou A = ±i 


Portanto a matriz é diagonalizável e pode-se tomar 




"00 " 




J = 


i 








-i 







Claramente, 



vi 



10 
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é um vector próprio de X = 0. Os vectores próprios de X = i são os que verificam 






a 




' " 




b 


= 







c 








portanto uma base dos vectores próprios para X = i é dada por 



v 2 = 



Como A é uma matriz real os vectores próprios associados a —i são os conjugados 
dos vectores próprios de X = i. Uma base é dada por 



V 3 = 



Portanto uma matriz de mudança de base é dada por 

S = 



1 








1 1 





— i í 



(e) Os valores próprios de A são as soluções de 



-X 

-1 





-1 

-A 
1 



= 



^^ -A á + A-A = 
^^ A = 

Os vectores próprios de X = são os que satisfazem a equação 





-1 

-1 o 
1 

6 = 

—a — c = 



-1 









a 




' " 




b 


= 







c 








A dimensão do espaço próprio c/e é 1 portanto há um único bloco de Jordan. Assim 

J = 



1 
1 




Uma base dos vectores próprios associados a X = é 



Vl 
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A segunda coluna da matriz S obtém-se resolvendo a equação 
(A - 0I)v 2 = vi 

0-10 

-1 -1 

1 o 





a 




' 1 




b 


= 







c 




-1 



Pode-se tomar 



V 2 = 



O 



-1 



O 



A terceira coluna da matriz S obtém-se resolvendo a equação 

(A - 0i> 3 = l>2 

0-10 

-1 -1 

1 





a 









b 


= 


-1 




c 








Pode-se tomar 



V:í 



Portanto uma matriz de mudança de base é dada por 

1 1 

s = 



-1 o 

o o 



D 



